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2 ■ Abstract 

-4— > . 

C^ , We calculate the leading-order matrix element for exclusive decays of 6 ^ S7 

Mh, in the quenched approximation of lattice QCD on a 24^ x 48 lattice at /3=6.2, 



using an 0(a)-improved fermion action. The matrix element is used to extract 



KA , the on-shell form factor Ti{q =Q) for B -^ K*^ and Bs -^ 07, using two 

H : 

^ different assumptions about the q^ dependence of the form factors for these 

decays. For B -^ K*-f, Ti{q'^=0) is determined to be 0.159l|^ ± 0.067 or 
0.124_i8 ± 0.022 in the two cases. We find the results to be consistent (in the 
Standard Model) with the CLEO experimental branching ratio of BR (B -^ 
K*j) = (4.5 ± 1.5 ± 0.9) X 10-5. 
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I. INTRODUCTION 

A. The Standard Model and New Physics 

Theoretical interest in the rare decay B — » K*'y as a test of the Standard Model has 
been renewed by the experimental results of the CLEO collaboration ||l[]. For the first time, 
this mode has been positively identified and a preliminary determination of its branching 
ratio given. 

The radiative decays of the B meson are remarkable for several reasons. The decay 
B — > K*'y arises from the fiavour- changing quark-level process b -^ 57, which occurs through 
penguin diagrams at one- loop in the Standard Model. As a result, the decay is a purely 
quantum effect and a subtle test of the Standard Model. The process is also sensitive to new 
physics appearing through virtual particles in the internal loops. Existing bounds on the 
b ^ S'j branching ratio have been used to place constraints on supersymmetry (SUSY) 0-|^ 
and other extensions of the Standard Model (SM) @,|10|- A comprehensive review of these 



results can be found in |Tl[]. Finally, it is also remarkable that this rare process has a 
sufficiently large branching ratio to be detected experimentally. Thus, accurate experimental 
measurements and accurate theoretical calculations of these decays could soon probe new 
physics at comparatively low energies. 

In order to compare the experimental branching ratio with a theoretical prediction it is 
necessary to know the relevant hadronic matrix elements. These have been estimated using 
a wide range of methods, including relativistic and nonrelativistic quark models 1T^-|I^ 



two-point and three-point QCD sum rules |15-20| and heavy quark symmetry |21], but 



there remains some disagreement between the different results. It is therefore of interest to 
perform a direct calculation of the matrix elements using lattice QCD. The viability of the 
lattice approach was first demonstrated by the work of Bernard, Hsieh and Soni |2^ in 1991. 



Excluding QCD contributions, the free quark decay 6 — > 37 in the SM proceeds by 
diagrams similar to that shown in Fig.(|I]). The charm and top quark dominate, because the 



up quark contribution to the loop is suppressed by the small CKM factor |V^feV^*s|. 

If the value of the top mass is assumed, the Standard Model can be tested by deriving 
an independent result for BR {B -^ K*'~f) . Deviations from the expected branching ratio 
would be an indication of contributions to the decay from physics beyond the SM, to which 
this decay is potentially sensitive. 

Research on such contributions can be classified into supersymmetric and non- 
supersymmetric extensions of the SM. In the latter case, Cho and Misiak ^^ considered 
SU{2)l ® SU{2)r left-right symmetric models and found considerable variations from the 
SM result for a wide range of the free parameters, while Randall and Sundrum p4| found 
significant potential deviations from the SM in technicolour models. Anomalous WW'y cou- 
plings in 6 — i> S7 have been analysed and the results found to be consistent with the SM. The 



bounds obtained from this approach can improve on those from direct searches |25-28|. The 



contributions from two Higgs doublet models p9|j30| have been analysed to obtain bounds 



on the charged Higgs mass and tan/3, the ratio of the vacuum expectation values of the 
doublets IPT1,|32 . 



SUSY models also involve additional Higgs doublets, but the contribution of other boson- 
fermion loops, in particular charginos (x^) with up type squarks, and gluinos (g) or neu- 
tralinos (x°) with down type squarks must also be included [3H8|J33||. A thorough study of 
the decay in the Minimal Supersymmetric Standard model can be found in reference [Q. 
There are a strong contributions from chargino and gluino loops, especially for large tan /3, 
which interfere destructively with the Higgs contribution and allow SUSY to mimic the SM 
in some regions of parameter space. As a result, the current limits on tan (3 and Higgs masses 
are weak, but will tighten as more stringent bounds on superpartner masses are obtained. 

For the rest of this paper, we shall use the SM as the appropriate model, and look for 
possible deviations from the experimental branching ratio. It should be noted that the 
lattice calculation is needed only to determine the effects of low energy QCD, and these are 
independent of new physics. The low energy effect of many extensions of the SM will be 
completely contained within the renormalisation group operator coefficients, and hence it is 



straightforward to allow for contributions from different models. 

B. Exclusive vs. Inclusive decay modes 

The inclusive decay B -^ Xg'j is predominantly a short distance process and can be 
treated perturbatively in the spectator approximation. It is also possible to use Heavy Quark 
Effective Theory (HQET) to compute the leading I /ml corrections ||3^. The experimental 



inclusive branching ratio has been determined at CLEO [35 



BR {B -^ X,7) = (2.32 ± 0.51 ± 0.29 ± 0.32) x 10"^ (1) 

The procedure for obtaining this result has a mild model dependence (the final result is a 
function of rrif,). 

In addition, the branching ratios of the exclusive decay modes oi b —>■ s'-f can also be 
experimentally determined, and the present published branching ratio for B —>■ K*'j from 
the CLEO collaboration [Q is, 

BR [B -^ K*-f) = (4.5 ± 1.5 ± 0.9) x 10"^ (2) 

The advantage of this measurement is that the experimental number is model independent. 
Theoretical calculation of the relevant exclusive matrix element requires the determination 
of long distance QCD contributions which cannot be determined perturbatively, but can be 
computed using lattice QCD. 

C. The Effective Hamiltonian and Hadronic Matrix Elements 

In order to determine the low energy QCD contributions to this decay, the high energy 
degrees of freedom must be integrated out, generating an effective AB = —1, AS = 1 



hamiltonian. Grinstein, Springer and Wise |^ determined the hamiltonian Ti-eff, to leading 
order in weak matrix elements. 



AC ^ 

neff = --jl-v,,v*Y.c,if^n, (3) 

where the eight operators Oj are multiphed by the renormahsation group coefficients, Ci{fi). 
Six of the operators are four-quark operators and two are magnetic moment operators, 
couphng to the gluon and photon |^. The operator which mediates the 6 — > 57 transition 



IS, 



<^7 = T^rn,sa,,hl + 75)6 F^^ (4) 

The coefficients Ci{^) are set by matching to the full theory at the scale /i = M^/. The coef- 
ficient Cj{mb) is determined using the renormalization group to run down to the appropriate 
physical scale ft = nib |^, and is approximately given by. 



where, in the Standard Model [39 



X /2 2 5 7 x(3x-2), \ mj 

- -, 7T7 \ -X H X ; — log X , X = 7i 

2X-13V3 12 12 2 a;-l ^ /' M}, 



CriMw) = T^TZ-^, io^' + TT^^-TT^- n^TZ-^^ogx] , x = -f . (6) 



The effects of scale uncertainty in the leading order approximation have been considered by 



Buras et al. [31 



To leading order, the on-shell matrix for B -^ K*'j is given by, 

M = '^0C,{m^)V,,VZv^*{KV,\B), (7) 

where, 

JfM = s(T^,uq''bR, (8) 

and rj and q are the polarization and momentum of the emitted photon. As outlined by 
Bernard, Hsieh and Soni ^2\, the matrix element {K*\scrf^,jq'^bii\B) can be parametrised by 



three form factors, 



{k*\j,\b) = J:c;tm'), (9) 



where, 



Cl = 2£^,,,eVfc'', (10) 

Cl = e^^iml -ml,)-e-q{p + k)^, (11) 



^3 _ . ^ L ^ 
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Cl = e.qlq,- , " , (p + A;) J ■ (12) 

As the photon emitted is on-shell, the form factors need to be evaluated at q^=0. In this 
hmit, 

T2(g2=0) = -tT,{q'=0), (13) 

and the coefficient of r3(g^=0) is zero. Hence, the branching ratio can be expressed in terms 
of a single form factor, for example, 

BR{B -. 1^7) = ^.rnlGlmlrs (l - ^] \V,,V:f\C,{m,)m{q'=0)\\ (14) 

This paper concerns the evaluation of Ti(0). We shall outline how matrix elements of 
the form {V\J^\P), where \P) is a heavy-light pseudoscalar meson and \V) is a strange- 
light vector meson, can be calculated in lattice QCD and explain the computational details 
involved. We shall evaluate the form factors Ti(g^=0) and T2(g^=0), make some statements 
about the systematic error, and compare the calculated value of BR [B -^ K*'y) with the 
results from CLEO. 

D. Heavy Quark Symmetry 

We cannot directly simulate 6-quarks on the lattice, as will be explained below. In- 
stead, we calculate with a selection of quark masses near the charm mass. This means that 
any results for the form factors must be extrapolated to the &-quark scale. Heavy quark 
symmetry PD[ tells us that. 



.2 N _, _l/2 



TlilmaxJ ■" ""P ,._. 

(15) 



imax) '^ '"P 



as the heavy quark mass, and hence the pseudoscalar meson mass, mp, grows infinitely large. 
Combining this with the relation T2{q^=0) = — iTi(g^=0) constrains the q'^ dependence of 
the form factors. 

Pole dominance ideas suggest that, 



(1 - qymjy 



^^(^^) - .. :z^.n. (16) 



for i = 1, 2, where rrii is a mass that is equal to mp plus 1/mp corrections and rii is a power. 
Since 1 — q^rnax!''^} ~ l/mp for large mp, the combination of heavy quark symmetry and 
the form factor relation at g^ = implies that rii = n2 + 1. Thus we could fit T2{q^) to a 
constant and Ti(g^) to a single pole form or fit T2{q^) to a single pole and Ti(g^) to a double 
pole. These two cases correspond to, 

{mp single pole 
3/2 • (1^) 

m.p double pole 

As we will see, our data for T2{q^) appear roughly constant in q^ when mp is around the 
charm scale, but have increasing dependence on q^ as the heavy quark mass increases. We 
will fit to both constant and single pole behaviours for T2(g^) below. 

II. LATTICE FIELD THEORY 

The hadronic matrix element iy\J^\P) for the b —^ s'-f transition can be obtained from 
the correlator {0\JJ^{x)T^^{y)Jp{0)\0), where Jp and J^ are interpolating fields for the P 
and V mesons, consisting of a heavy quark, h, a light quark, /, and a strange quark, s; 

Jp{x) =l{x)^5h{x), (18) 

J^{x)=J{x)-fps{x), (19) 

T,,iy) = siy)a^My)- (20) 

The full matrix element (V|so"^,y|(l + ^rj)h\P) can be derived using the Minkowski space 
relation, 
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f^y'" = 2^'^'''^PA- (21) 

In Euclidean space, the correlator (0| J^(x)T^j,(|/) Jp(0)|0) can be computed numerically 
using the functional integral, 

(0|j;(x)r^,(l/)4(0)|0) = ^jVAVqVq j;(x)T^,(y)4(0) exp(-5[A, g, g]), (22) 

= ^j'DA Tr(75ff(0,i/)a^,5(i/,x)7,L(x,0))exp(-^eif), (23) 

where S[A, q, q] is the QCD action and S{y, x), H{y, x), L{y, x) are the propagators from x 
to y for the s, h and / quarks. Working in momentum space, we calculate the three-point 
correlator, 

Cf,i{t, tf, p, q) = Y: e''^-'e-^'^y{Jp{tf, x)T,,{t, f)4,(0)) (24) 



x,y 



t,tf-t,l ^oo ^ ZrLp ZrLy 

To obtain the matrix element {P{p)\hafj_i^s\V (k)) , we take the ratio, 

Cp {tf-t,p)Cy'' {t,p-q) 
where the two-point correlators are defined as, 

X 



2E 



(e-'^pt ^ ^-E,iT-t)^^ ^ (27) 

Cr{t.k) = -(^) Ee^'-'(4.(i,^)-/y(0)) 

2Ev ^ ' 



By time reversal invariance and assuming the three points in the correlators of Eq. (p6|) are 
sufficiently separated in time, a term proportional to the required matrix element dominates: 

Cpp.^ —? ^=^Y.^p{V{k,t)\-sa^,h\P{p))^..., (29) 

t,tf~t,l — >oo ZjpZjy 



and Cp^y approaches a plateau. The three-point correlator is calculated in its time reversed 
form to allow the use of previously calculated light propagators. The factors Zp, Zy and 
the energies of the pseudoscalar and vector particles are obtained from fits to the two-point 
Euclidean correlators. 

In order to simulate this decay on a sufficiently finely spaced lattice, vacuum polarisation 
effects were discarded and the gauge field configurations generated in the quenched approx- 
imation. The decays D -^ Keu, D -^ K*ev, and Ds — > 0ez/ have been calculated in the 



quenched approximation [p| -^ and have been found to be in relatively good agreement 
with experiment. It is therefore quite plausible to assume that the systematic error from 
the quenched approximation for this calculation would be of a similar size. 

The matrix element {K*\'safj^,jq'^h^\B) cannot be directly calculated as realistic light 
quarks cannot be simulated owing to critical slowing down in determining the propagator for 
small masses. Instead light quarks are simulated at a number of masses approximately that 
of the strange quark mass, and any result is extrapolated to the chiral limit. Furthermore, 
h quarks cannot be simulated directly as the 6-quark mass is greater than the inverse lattice 
spacing (2.73(5) GeV), and the variation of the propagator would occur over lengths smaller 
than the lattice spacing. As a result, heavy quarks are simulated with masses around the 
charm quark mass, and the results extrapolated to m;,. Hence, {V\Jfj\P) has to be calculated 
at a number of different light, strange and heavy quark masses. 

III. COMPUTATIONAL DETAILS 

Sixty SU{?>) gauge configurations were generated in the quenched approximation for a 
24^ X 48 lattice aX (3 = 6.2. These configurations were generated with periodic boundary 
conditions using the hybrid over-relaxed algorithm, and the standard discretised gluon ac- 
tion, defined in |^^. The configurations were separated by 400 compound sweeps, starting 
at sweep number 2800. The inverse lattice spacing was determined to be 2.73(5) GeV, by 
evaluating the string tension [H3. In physical units, this corresponds to a spacing of ap- 
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proximately 0.07 fin and a spatial size of 1.68 fm. In order to simulate heavy quarks whose 
masses are approaching the inverse lattice spacing, the 0(a)-improved fermion action of 
Sheikholeslami and Wohlert [^ (also referred to as the clover action) was used. This is 
defined as, 

Sp = S^ -i-Y^ q{x)Ff,^{x)a^^q{x), (30) 

where S^ is the standard Wilson fermion action [^,0 and -F^,y is a lattice definition of the 



field strength tensor, which we take to be the sum of the four untraced plaquettes in the fj,^ 
plane open at the point x, 



F^u{x) = 7 I] TT Unf.u{x) - f/^Jx) 



4 ^1 2t 



n/i!/\ 



(31) 



In using this action, all observables with fermion fields q, q must be "rotated" , 



q{x) -. (l-i^)g(x), (32) 

q{x) -^ q{x)(^l + -^y 

where A^ is the discretised covariant derivative, operating on the quark fields as, 

^t^ <l{x) = -(Uf,{x)q{x + fi) - Ul{x- fi)q{x- fi)y (33) 

q{x) A^ = - (q{x + fi)Ul{x) - q{x - fi)U^{x - /i)) . 



This action eliminates the tree level 0(ma)-error of the Wilson action [^], which can be 
significant for heavy quark systems [^9| , ^ . 

For each configuration, quark propagators were calculated using the over-relaxed minimal 
residual algorithm with red-black preconditioning for k, = 0.14144, 0.14226 and 0.14262, 
using periodic boundary conditions in the spatial directions and anti-periodic boundary 
conditions in the temporal direction. Smearing was not used in the calculation of these light 
propagators. The first two n values can be used to interpolate to the strange quark mass 
which corresponds to k = 0.1419(1) [^. The third n value, corresponding to a somewhat 
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lighter quark, was used in conjunction with the others in order to test the behaviour of the 
data in the chiral hmit. 

Heavy propagators, for Kh = 0.121, 0.125, 0.129 and 0.133, were evaluated using timeslice 
24 of some of the above propagators as the source. For Kh = 0.121 and 0.129, the propagators 
for all of the light n values were used. For k^ = 0.125 and 0.133, the propagators for 
K = 0.14144 and 0.14226 were used. To reduce excited state contamination, these sources 
were smeared using the gauge invariant Jacobi algorithm [^, with an r.m.s. smearing radius 



of 5.2. Because of memory limitations, these propagators were evaluated only for timeslices 
7 to 16 and 32 to 41. 

Using these propagators, the three point correlators were evaluated. The spatial momen- 
tum p was chosen to be (0, 0, 0) or (vr/12, 0, 0) (the lowest unit of momentum in lattice units 
that can be injected). All possible choices of q were calculated such that the magnitude of 
the spatial momentum of the vector meson k was less than v27r/12. This is because the 



signal of light hadrons degrades rapidly as the momentum is increased p3 |. 

In order to obtain (VlscXf^uhlP), the decay constant and energy were determined for the 
pseudoscalar of each heavy-light k combination and the vector of each possible light k, com- 
bination, for all possible momenta used. The process of extracting these is well understood 
and has been discussed in detail elsewhere [^. As the two point functions are periodic, a 



correlator at a time < t < 24 was averaged with the same correlator at 48 — t to improve 
the statistical sample. This "folded" data was fitted to Eg . (|27D or Eq.(^) for timeslices 15 
to 23. For both the two point and three point functions we utilised the discrete symmetries 
C, P and T (folding) wherever possible, in addition to averaging over equivalent momenta. 



The statistical errors for all correlators were determined by the bootstrap procedure [gj] , us- 
ing 1000 bootstrap subsamples from the original configurations. The finite renormalization 
needed for the lattice-continuum matching of the a^u operator has been calculated |^ but 
has a negligible effect here (0(2%)) and was not included. It introduces a small correction 
to the branching ratio which is considered in the conclusions. 

As outlined in the previous section, the weak matrix elements Cpf^i, were extracted from 

11 



the three point data and the fits to the two point data. Having divided out the contributions 
from the two point amphtudes and energies, the matrix element {V\safj,^h\P) was isolated. 
These matrix elements were combined to determine the form factors Ti(g^), T2(g^^^) and 
T2(g^). Each form factor was extracted by a correlated fit to a constant for timeslices 11, 12 
and 13. 

IV. RESULTS 

The data for unphysical masses, and off-shell photons must be combined to isolate the 
form factors and extrapolate to the physical regime. It is clear from Eq.(0) and Eq.(|l4D 
that the branching ratio can be evaluated from Ti(g^=0; m^; rrii^.) or T2(g^=0; m^; rrii^.). 
As demonstrated in a previous paper |^, the evaluation of Ti(g^=0; mp; rriK*) is relatively 
straightforward, and T2{q'^=0;'mp;mK*) can be determined in a similar way. To test heavy 
quark scaling, we also extracted the form factor T2 at maximum recoil, where q^ = q^ax — 
{mp — mvY, in the same way as Bernard et al. |^^. These form factors were extrapolated 
to the physical mass mp = m^B, and an estimate of systematic errors in the extrapolation 
made by comparing different methods. 

A. Extraction of form factors 

1. ri(g2) 

The form factor Ti can be conveniently extracted from the matrix elements by considering 
different components of the relation, 

4(fcV - Kfc^)Ti(g') = e'^^'^Cp^^q". (34) 

We see a plateau in Ti about t = 12. The use of smeared operators for the heavy quarks 
provides a very clean signal, with stable plateaus forming before timeslice 11. The data for 
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the heaviest of our hght quarks, ni = Kg = 0.14144, with the smallest statistical errors, are 
shown in Fig.(^. 

The form factor is evaluated for each of the five possible values of g^. We fit Ti(g^) to a 
pole or dipole model in order to obtain the on-shell form factor Ti(g^=0), 

r.(,^).^(^. r.tf ) ^ -^(4z2) . (35) 

1 — q^/m^ (1 — q-^ /m'^Y 

We allow for correlations between the energies of the vector and pseudoscalar particles and 
Ti at each q^. An example of such a fit, for ki = Kg = 0.14144, is shown in Fig.(^ and the 
full set of fit parameters and their x^/d.o.f. are shown in tables |I| and 0. 

The chiral limit behaviour of Ti(g^=0; mp; mv), interpolated from a single pole fit, was 
explored for Hh = 0.121 and 0.129, in our earlier work |]56|. To test for approximate spectator 
quark independence, we compared the single pole fits of the form factor to the two functions, 

Ti{q'^=0;mgjight) =a + brrii, (36) 

Ti{q'^=0;mgjight) = c, (37) 
where mi is the lattice pole mass, 

mi = l-{--^), (38) 

and Kcrit = 0.14315(2) |5^. The linear coefficient b was found to be consistent with zero 
for each combination of Hg and Hh (see Fig.^. From table 0, the x^/d-o.f. for both fits 
are similar, indicating that for the data available, the assumption that the form factor is 
a constant, independent of the spectator quark mass, is valid. Hence, the data for ni = 
0.14144 was used for the chiral limit, and a simple linear interpolation carried out between 
Kg = 0.14144 and 0.14226 for the strange quark, in order to obtain Ti{q'^=0;mp;mK*)- 
These results are listed in the columns labelled (b) and (c) in table |V114 

2. Uq^) 

The form factor T2 can be extracted from the matrix elements using the same procedure 
as Ti, by considering the different components of, 
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{nip - my)T2{q^] mp] my) = Cu^q", (39) 

for all i (not summed) such that g* = 0. A typical plateau for T2 is shown in Fig.(^. We 
extract T2 for a range of g^ as shown in Fig.(^. 

Fig.(l^) shows that T'2(g^) is roughly constant as a function of q^ for our data, with heavy 
quark masses around the charm mass. We fit T2 to a constant: we can then compare with 
the value of Ti(g^ = 0) where Ti is fitted with a single pole form. We also fit T2 to a single 
pole form (as shown in the figure) and compare with Ti(g^ = 0) when Ti is fitted with a 
double pole form. The results of the fits for T2 are shown in tables |T| and |I^, and the chiral 



extrapolations for the single pole fit in table 0. The pole mass is found to be large, and 
a linear behaviour holds well for all possible g^, including qf^ax^ ^s shown in Fig.(|^). Once 
again the data for ki = 0.14144 was used for the chiral limit and the results are listed in the 
columns labelled (d) and (e) in table |V111| . 

The ratio Ti(g^=0; mp; mi^.)/T2(g^=0; mp; m^.) is shown in Fig.(|^). The two sets of 
points show Ti fitted to a double pole form and T2 to a single pole or Ti fitted to a single 
pole and T2 constant. The ratio should be 1, in accordance with the identity T2(0) = 
— iTi(O), Eq. (|l3|) . We find greater consistency from the double-pole/single-pole fit. 



3. T2[q, 



max) 



The evaluation of T2(g^Q^; mp; my) is also straightforward, since at zero momentum, 
p=0, k=0, the contributions from other form factors vanish, 

[mp + my) T2{ql^ax) = '^110 (p = 0, k = 0), 

= C220(p = 0,k = 0), 

= C33o(p = 0,k = 0). (40) 

An example of this data is shown in Fig.(^). The behaviour oiT2{q^ax'i '^^p'i ^^v) as a function 
of the spectator quark mass was examined at ^h = 0.121 and 0.129 in the same way as for 
Ti(g^=0). It was again found that the linear coefficient b was consistent with zero for each 
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combination of Kg and Hh'- see Fig.(P) for an example. From table |VI1| , the x^/d.o.f. for 
both fits are seen to be similar, indicating that for the data available, the assumption that 
the form factor is independent of the spectator quark mass is valid. Hence, the data for 
Ki = 0.14144 was used for the chiral limit, to obtain T2(g^„3..; mp; rriK*)- 

Bernard et al. |E^ converted this result to q^=0 by assuming single pole dominance. 



The current J^ in the matrix element can be expressed in & V + A form, with Ti corre- 
sponding to the vector component and T2 and T3 to the axial current. Therefore, in a 
single pole model, the exchanged particle, Pgi, for the T2 form factor should be the low- 
est J^ = 1"^ state with the correct spin, parity and strangeness quantum numbers. We 
extracted T^'"^(g^=0; mp; m^^.) from T2{q^ax) using a single pole model, with the mass of 
the 1^ states determined from fits to two-point functions for each heavy quark mass. The 
results of these extrapolations are shown in the column labelled (a) in table |V111| . 

The ratio Ti{q'^=0; mp; mK-) /T2° '^{q'^=0; mp; rriK*) is shown in Fig. (p!0|) . We note that 
using a fixed pole mass from two-point functions gives a 10-20% difference in the ratio (at 
the heaviest masses) compared with allowing the pole mass to vary in the fits. 

B. Extrapolation to Mb 

The appropriate ansatz for extrapolating the on-shell form factor in the heavy quark mass 
to Ti(g^=0; rriB', rriK*) is not a priori clear. As we saw in section [1D| , one has to model the 



q^ dependence of the form factors, maintaining consistency with known heavy quark scaling 
results [|^ at q^ax^ from Eq. ([T5|) , and the relation Ti(0) = iT2{0). Expanding unknown 



parameters in powers of l/rrip, one obtains scaling laws for the on-shell form factors Ti(g^=0) 
and T2(g^=0). Thus, while the scaling behaviour of T2(g.^„^) can be checked directly, the 
behaviours of Ti(0) and T2(0) will depend on assumptions made for the g^ dependence. We 
now address these issues. 
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Bernard et al. |5^ used the heavy-quark scahng law for the off-shell form factor, 
T2{q'^^^;mp;mK*) to extrapolate T2 to T2(g^^^; m^; mi^-.), before applying a single pole 
dominance model as before to reach the on-shell point T2(g^=0; ms', uik*)- They estimated 
the appropriate pole mass. The validity of the pole model over the wide range of momentum 
transfer from g^=0 to q^^x ^^^ required, but tests at heavy quark masses around the charm 
quark mass showed it to be quite accurate. 

Our results for T2(g^; mp] rriK*), see Fig.(0), appear nearly independent of q^ for masses 
mp around the charm scale. Hence, we have fitted T2 to both single pole and constant forms, 
with corresponding behaviour for Ti. This will give us two alternative forms for the heavy 
mass dependence of Ti(g^=0; m^; niK*). 



1- T2{q, 



2 ) 

max/ 



At q =qmaxy ^^^ initial and final hadronic states have zero spatial momentum and the 
contributions of form factors other than T2 vanish, 

{K*\saf,^q''bR\B) = e^.(m| - m^.)T2(g^„^). (42) 

In the heavy quark limit, the matrix element of Eq.p^) scales as m^ , owing to the normal- 



isation of the heavy quark state [y/rriB) and the momentum transfer q (g° = mp — mx*)- 
The leading term in the heavy quark scaling of T2(g^„^.) is expected to be m^ ' , analogous 
to the scaling of fp Il59|j49(| . Higher order 1/mp and l/rn^p corrections will also be present, 
as will radiative corrections [|60| , |6l|] . 

Hence, the form factor T2{q^ax) should scale as. 



T2{qmax^ TTip] mK*)Vrri^ = const, x [as{mp)] ^/^^ 1 H H 1- + 

Y mp rrip 

To test heavy quark scaling, we form the quantity, 

- , 2 \ /^p f as{mp)\ 

y mp \as{mp)J 

where we approximate ^^(/i) by, 
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(43) 



(44) 



with Aqcd = 200 MeV and Pq = 11 — |A^/- In the quenched approximation, Nf is taken 
to be zero. The normahsation ensures that T2 = ^2(9^^3.) at the physical mass ms- Linear 
and quadratic correlated fits to Eq.(^) were carried out with the functions, 

f2{mp) = A(l + —), (46) 

V nip J 

f,(»,,=^(l + ^ + ^), (47) 

and are shown in Fig.(|TT]). Taking the quadratic fit of T2 at mp = tub as the best estimate, 
and the difference between the central values of the linear and quadratic fits as an estimate 
of the sytematic error, T2 was found to be 

T2{ql,ax\ ^B\ mx*) = 0.269+^^ ± 0.011. (48) 

Once 72(g^„^) is extracted, we can obtain T2(0) in the two cases, pole model or constant, 
for the (^ behaviour. If T2 is constant, then Eq. (^8D is the result at q^ = 0. In the pole 
model, the expected exchange particle for T2 is the 1+ -B^i state, but experimental data for 
its mass is not yet available. However, it is possible to estimate reasonable upper and lower 



bounds for the mass from HQET. It can be shown that [59| 



— A 1 

mB,,-mB = AA + — + 0(^), (49) 

nib ml 

— A 1 

mD.i -m^ = AA + — + 0(^). (50) 

nic m^ 

Neglecting terms of order l/ni^, the upper and lower bounds for Eq. (|49|) are, 

—^{rriDsi - rno) < niB,-, - niB < nio,^ - nio (51) 



Making the approximation, 

rric mp, + Smu* 



^6 miB + Sm^* 
the range of the expected pole mass can be found, 
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(52) 



ms^^ = 5.74 ±0.21 GeV. (53) 

Therefore, 

Tr'\q'=Q-, mj,- ruK*) = O.lUtUlt, (54) 

where the first error is statistical and the second is the systematic error obtained by combin- 
ing the variation of the pole mass within its bounds and the systematic error from Eq.(^). 
There is clearly a significant systematic difference between the results in Eq.(^8D and Eq.(p4D 
corresponding to the two assumed forms for T2(g^). 

2. ri((?2=o) 

If constant-in-g^ behaviour is assumed for T2, then T2(0) should satisfy the same scaling 
law as T2(g^„^) in Eq. (|4^) . Combining this with the identity T'i(O) = iT2{0) leads to a 
scaling law for T'i(O): 



Ti(0: mp: mK*)x/rn^ = const, x [a Jmp)l"^/^° (l + — + -^ + ...). 

y mp nip J 



(55) 



If single pole dominance is assumed for T2 and the mass of the exchanged 1+ particle 
can be expanded as, 

mp^^ = mp IH — H \- + . . . , (56) 

y mp rrip J 

then Ti(g^=0; mp; m.^*) should satisfy a modified scaling law, 

Ti(0; mp- mK*) rn/^ = const, x [a,(mp)]-2/^° ( 1 + — + ^ + . . . ) , (57) 

Y m,p irip ) 

where the unknown coefficients in Eq.(^) have been absorbed into the unknown scaling 
coefficients of the matrix element. A similar scaling relationship has been found by Ali et 
al. [|1^] by the sum rules approach. 

The two scaling forms were tested in the same way as for T2(g^„^), by forming the 
quantities. 



f,.r,„M,(-)"'^(-|-))^'r (58) 

KniB/ \as{mB)J 

where A^ is 1 or 3 as appropriate. 

Linear and quadratic fits were carried out with the same functions as for T2, allowing for 
correlations between masses and form factors. They are shown in Fig. ([12|) . The x^/d.o.f. 

'5/2 

was approximately 1 for the rrip scaling law, indicating that the model is statistically valid 
in the available mass range. For the nip scaling law we found a x^/d.o.f. of 0.3. 
The correlated quadratic fit with radiative corrections gives, 

(0.l59tlt m}/%caling 
Ti{q =0]mB]mK*) = { , (59) 

[ O.UAtll nif scaling 

where the errors quoted are statistical. 

All methods of evaluating Ti{q'^=0;mp]mK*) at intermediate masses are compared in 



table [VII1| . We consider the differences between the methods as a measure of part of the 
systematic error. The differences between the methods of determining the form factors at 
the computed masses are of a similar size (~ 10%) to the systematic error at the physical B 
mass, as measured by the linear or quadratic extrapolation of Ti in the inverse heavy meson 
mass. 

The final result for Ti(g^=0; ttib', rriK*) is taken from the quadratic fit for Ti, with an esti- 
mated systematic error in extrapolation given by the difference between linear and quadratic 
fits, 

[ 0.1591^^ ± 0.067 m](^ scaling 
Ti{q^=0;mB;mK*) = \ ■ (60) 

[ 0.1241?^ ± 0.022 m^/^ scaling 

The extrapolation is shown in Fig.(0). We note that the value obtained from rrip scaling is 
consistent with the corresponding value from T2 calculated using the single pole g^ behaviour 
discussed earlier. 
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C. Bs -^ (t>-f 

Much of the analysis above can also be applied to the rare decay Bg -^ 07. ALEPH 
and DELPHI |63| have looked for this decay and obtained 90% CL upper bounds on its 
branching ratio of 4.1 x 10~^ and 1.9 x 10~^ respectively. Future research into this decay 
at LEP is planned. The branching ratio for this decay can be expressed in a form similar 
to Eq.(|l^ 



C2 \ 3 
l-£^j \V,,V:f\C,{mt)miq'=0)\', (61) 

where Tf is the relevant form factor from the decomposition of (0| J^|i?s). In determining 
this matrix element numerically, the interpolating operator Jj (x) is replaced by the operator 
Jp{x) defined as, 

J}{x) = s{xhps{x). (62) 

As a result of the presence of two identical particles in the final state, there is an extra 
additive term in the trace of Eg. (p3D , which corresponds to ss creation from purely gluonic 
states. It is expected that this process is heavily suppressed by Zweig's rule [|6^ -|66[|, and 
hence the extra term is neglected. 

As the variation of the form factors with respect to the spectator quark mass has been 
discarded, it can be assumed that, 

Tl{q^=0]mp;m^) = Ti(g^=0; mp; mi^-J, (63) 

^2(5^=0; mp] m^) = T2(g^=0; mp] mx*)- (64) 

By employing the same ansdtze for extrapolating Ti and T2 as the previous sections. 

f 0.165^30 ± 0.060 m!/^ scaling 
mq'=0;ms^;m,) = { "'^ , (65) 

[ 0.125t?^ ± 0.021 mp/^ scaling 
T^{qL., mB/,m^) = 0.2701^7 ± 0.009, (66) 

Tr%'=0; ms^ m^) = O.lUtltH. (67) 

We note that Tl(q'^=0), with rrip scaling, and T2*''^°'^(g^=0) are consistent with each other. 
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V. CONCLUSIONS 

In this paper we have reported on an ab initio computation of the form factor for the 
decay B -^ K*'~f. The large number of gauge configurations used in this calculation enables 
an extrapolation to the appropriate masses to be made and gives a statistically meaningful 
result. To compare this result with experiment we convert the preliminary branching ratio 
from CLEO, BR{B -^ K*-f) = (4.5 ± 1.5 ± 0.9) x lO^^ based on 13 events H, into its 
corresponding Ti form factor, assuming the Standard Model. We work at the scale fi = mi, = 



4.39 GeV, in the MS scheme, using a pole mass of Mi, = 4.95(15) GeV [q^I to determine 
mb ||6^. Taking iV^s^tbl = 0.037(3) |6^, tb = 1.5(2) ps |70| , |7T| and all other values from the 



Particle Data Book combined with Eq.([T|), we find T^''^ to be 0.23(6), 0.21(5) and 0.19(5) 
for top quark masses of rrit = 100, 150 and 200 GeV respectively. We find the calculated 
value for Ti consistent with these results to within two standard deviations. 

In calculating the branching ratio, we use the perturbative renormalisation of a^^, ||5^ 
with a boosted coupling, g"^ = l.Tf^Q, and the anomalous dimension, 'jg^q = — (8/3)((7^/167r^), 
to match the lattice results to the continuum at the scale ^ = mi,-, giving a matching 
coefficient oi Z k, 0.95. We apply a correction of Z^ = 0.90 in the calculations below. 
Varying the scale of C-j{^) from /i = ■025/2 to yU = 2mi, changes the final branching ratio 
by +27% and —20% respectively. This is due to the perturbative calculation of Cj{^) 
and future work on next-to-leading logarithmic order corrections will reduce this variation 
significantly ||3ll . 



These uncertainties cancel in the dimensionless hadronisation ratio, i?, 

which we find to be, 

' (l4.5l^^ (stat.) ± 6.1 (sys.) ± 1.6 (exp.) ) % m)/^ scaling 
< 
(8.8+i (stat.) ± 3.0 (sys.) ± 1.0 (exp.) ) % m^J.'^ scaling 
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R 



(70) 



Assuming the recent tentative result for rrit from CDF [^, the lattice results give a 
branching ratio for the decay B -^ K*'~f of, 

^2.5tn (stat.) ± 2.1 (sys.) ± 0.6 (exp.) tl (scale)) x 10"^ m)/^ scaling 
BR{B ^ K*-f) = { ) / ^ 

1.5tl (stat.) ± 0.5 (sys.) ± 0.3 (exp.) tt (scale)) x 10"^ m^/^ scaling 



(71) 

where we separate the statistical and systematic errors from the lattice, experimental and 
theoretical (scale) uncertainties. Combining errors to produce an overall result yields, 

(2.5 ± 1.3 (stat.) 1| (sys.)) x 10"^ mj/^ scaling 
BR {B -^ K*-f) = {) / . (72) 

1.5 ± 0.6 (stat.) tl (sys.)) x 10"^ m^f scaling 



we 



(74) 



Similarly for Bs -^ </)7, using m^^ = 5.3833(5) GeV [[f3| , [f^ and tb^ = 1.54(15) ps ||75 
find, 

^2.8t\l (stat.) ± 2.1 (sys.) ± 0.5 (exp.) tl (scale)) x 10"^ mj/^ scaling 
BR (Bs ^(p^) = { ) / , , (73) 

(l.6l^ (stat.) ± 0.6 (sys.) ± 0.3 (exp.) tt (scale)) x 10"^ rnf scaling 

(2.8 ± 1.2 (stat.) til (sys.)) X 10"^ m]^^ scaling 

(1.6 ± 0.6 (stat.) ^9° (sys.)) x 10"^ m^/^ scaling 

In obtaining these results, we have made some assumptions. Since this calculation is car- 
ried out with one lattice spacing, we cannot explore discretisation errors. However, the use of 
an 0(a)-improved action is expected to reduce these substantially. As the form factors and 
mass ratios evaluated are dimensionless, we also expect some of the systematic error from 
setting the scale to cancel. The extrapolation of matrix elements to the chiral limit has been 
neglected, although the current data indicates a weak dependence on the spectator quark 
mass. Without doing a simulation using dynamical fermions, the error due to quenching 
cannot be accurately estimated. However, the good agreement with experiment for other 
semileptonic, pseudoscalar to vector meson decays [^,^], that have been determined using 



coarser lattices and lower statistics, suggests that these errors are small. We find our re- 
sults consistent with previous calculations [^,0]. With form factors available over a range 
of masses, we have been able to incorporate heavy-quark symmetry into our extrapolation 
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and investigate phenomenologically motivated pole- dominance models. These methods su- 
percede the simple linear extrapolation used as a guide in our earlier preliminary study, 
where the limited set of two masses precluded an investigation of different extrapolation 
methods |5^. 

Whether pole dominance is a valid model for a large range of q^ is an important question. 
We have quoted results for two different possibilities for the q^ dependence of the form 
factors. Although the lattice results visually favour T2 constant in g^, at least for heavy 
quark masses around the charm mass, our fits favour a single pole vector dominance form 
for T2. The difference between the results indicates the need for a better understanding of 
the combined q^ and heavy quark scaling behaviour of the relevant form factors. 

We have not applied the constraint T'i(O) = 'iT2(0) to our fits in this paper, using instead 
the consistency of our results with this relation as a guide to the fitting method. We find 
that the the single pole dominance model for the q^ behaviour of T2 (and corresponding 
dipole behaviour for Ti) gives the most consistent fit. In this case we have attempted to 
determine the systematic consistency by comparing Ti(g^=0; tti^; mx*), extracted using the 
nip scaling law, with T2(g^=0; m^; ?7i^.) assuming pole model behaviour for T2 and the 
expected pole mass. It could be argued that both methods are equivalent. However, in 
extrapolating the form factor Ti{q'^=0;mp;mK*) to rriB, the coefficients in the fit are not 
fixed, which is equivalent to letting the pole mass vary. We require only that the leading 
order behaviour of Ti satisfy the mj dependence. 

We look forward to improved experimental results for the decay B — > K*'~f and observa- 
tion of Bs — > 07. We hope future lattice studies will significantly increase the accuracy of 
these calculations. 
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TABLES 



l^h 


Ks 


K/ 


low (ga)2 


high {qaf 


Ti(0) 


xVd.o.f. 


0.12100 


0.14144 


0.14144 


-0.032l^ 


0.258+^ 


0.283t}9 


11.2/3 


0.12100 


0.14144 


0.14226 


-0.035l^ 


0.265ti^ 


0.282t^; 


8.8/3 


0.12100 


0.14144 


0.14262 


-0.035lJ^ 


0.27ltl^ 


0.2941^^ 


4.1/3 


0.12100 


0.14226 


0.14144 


-QMA+f 


0.290+1^ 


0.27ll?^ 


7.3/3 


0.12100 


0.14226 


0.14226 


-0.0221^^ 


0.298l|^ 


0.279+^0 


6.7/3 


0.12100 


0.14226 


0.14262 


-0.028+1 


0.306+^^ 


0.294tn 


4.2/3 


0.12500 


0.14144 


0.14144 


-0.118l^ 


0.157t^ 


0.307^9^ 


10.3/3 


0.12500 


0.14226 


0.14144 


-0.104l^ 


0.183+5*0 


0.292+^^ 


8.5/3 


0.12900 


0.14144 


0.14144 


-0.188l^ 


0.086t^ 


0.335+^^ 


9.4/3 


0.12900 


0.14144 


0.14226 


-0.190l^ 


0.084t^ 


0.333li^ 


7.2/3 


0.12900 


0.14144 


0.14262 


-0.190l^ 


0.084t^ 


0.336+^? 


3.2/3 


0.12900 


0.14226 


0.14144 


-0.177l^ 


0.097+? 

— o 


0.3181^^ 


7.1/3 


0.12900 


0.14226 


0.14226 


-0.182l^° 


0.096^^2 


0.319+^^ 


7.5/3 


0.12900 


0.14226 


0.14262 


-0.186ti| 


0.103ti? 


0.324^37 


3.6/3 


0.13300 


0.14144 


0.14144 


-0.240l^ 


0.034t^ 


0.366l^^ 


7.3/3 


0.13300 


0.14226 


0.14144 


-0.233l^ 


0.041+^ 


0.343+^J 


5.7/3 



TABLE L Results of pole fits to Ti(q ; mp; my)- 
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Kh 


Ks 


Kl 


low (ga)^ 


high {qaf 


Ti(0) 


xVd.o.f. 


0.12100 


0.14144 


0.14144 


-0.032l^ 


0.258l^ 


0.279t?2 


9.2/3 


0.12100 


0.14144 


0.14226 


-0.035l^ 


0.265+1^ 


0.275+^^ 


7.8/3 


0.12100 


0.14144 


0.14262 


-0.035lJ^ 


0.27ltl^ 


0.2861^6 


3.5/3 


0.12100 


0.14226 


0.14144 


-0.014+^0 


0.290li^ 


0.266li^ 


6.4/3 


0.12100 


0.14226 


0.14226 


-0.022^^5 


0.298l|^ 


0.274+^^ 


6.5/3 


0.12100 


0.14226 


0.14262 


-0.028ii 


0.306ti6 


0.286t^^ 


3.9/3 


0.12500 


0.14144 


0.14144 


-0.118l^ 


0.157+^ 


0.308+9^ 


8.6/3 


0.12500 


0.14226 


0.14144 


-0.1041? 


0.183+5'o 


0.291+f^ 


7.3/3 


0.12900 


0.14144 


0.14144 


-0.188^4 


0.086t^ 


0.337l^^ 


7.9/3 


0.12900 


0.14144 


0.14226 


-0.190l^ 


0.084+^ 


0.334ti 


6.5/3 


0.12900 


0.14144 


0.14262 


-0.190l^ 


0.084+^ 


0.3371^2 


2.7/3 


0.12900 


0.14226 


0.14144 


-0.177l^ 


0.097+^ 


0.319+ig 


5.8/3 


0.12900 


0.14226 


0.14226 


-0.182+^° 


0.096ti2 


0.320li^ 


6.8/3 


0.12900 


0.14226 


0.14262 


-0.18611^ 


0.103+1^ 


0.325+^7 


3.2/3 


0.13300 


0.14144 


0.14144 


-0.240l^ 


0.034+^ 


0.362+^^ 


6.2/3 


0.13300 


0.14226 


0.14144 


-0.233l^ 


0.04ll^ 


0.34ll^° 


4.4/3 



TABLE II. Results of dipole fits to Ti{q'^;mp; my)- 
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Kh 


Ks 


Kl 


low (ga)^ 


high {qaf 


^2(0) 


xVd.o.f. 


0.12100 


0.14144 


0.14144 


-0.032l^ 


0.289l^ 


0.30ll?i 


6.4/4 


0.12100 


0.14144 


0.14226 


-0.035l^ 


0.293+5*0 


0.310+^^ 


7.5/4 


0.12100 


0.14144 


0.14262 


-0.035lJ^ 


0.297+1^ 


0.3151^4 


8.8/4 


0.12100 


0.14226 


0.14144 


-0.014+^0 


0.318+5'i 


0.288tlj 


4.4/4 


0.12100 


0.14226 


0.14226 


-0.022^^5 


0.324l|^ 


o.sootll 


5.6/4 


0.12100 


0.14226 


0.14262 


-0.028^1 


0.330ti 


0.312t^? 


5.6/4 


0.12500 


0.14144 


0.14144 


-0.118l^ 


0.190+^ 


0.322+19 


4.3/4 


0.12500 


0.14226 


0.14144 


-0.1041? 


0.214t^ 


0.310t| 


3.7/4 


0.12900 


0.14144 


0.14144 


-0.188^4 


o.iosti 


0.35ll^^ 


7.9/4 


0.12900 


0.14144 


0.14226 


-0.190l^ 


0.109+^ 


0.353t| 


6.9/4 


0.12900 


0.14144 


0.14262 


-0.190l^ 


0.112+^ 


0.3491^0 


10.7/4 


0.12900 


0.14226 


0.14144 


-0.177l^ 


0.126+? 


0.3321^^ 


5.2/4 


0.12900 


0.14226 


0.14226 


-0.182+^° 


0.129+?o 


0.337^20 


6.8/4 


0.12900 


0.14226 


0.14262 


-0.18611^ 


0.133+15 


0.336t^^ 


12.4/4 


0.13300 


0.14144 


0.14144 


-0.240l^ 


0.045+3 


0.372+^^ 


7.7/4 


0.13300 


0.14226 


0.14144 


-0.233l^ 


0.057l^ 


0.35ll^^ 


7.5/4 



TABLE III. Results of pole fits to T2{q ,mp; mv)- 
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Kh 


Ks 


Kl 


low (ga)^ 


high {qaf 


^2(0) 


xVd.o.f. 


0.12100 


0.14144 


0.14144 


-0.032l^ 


0.289l^ 


0.344t^^ 


18.4/5 


0.12100 


0.14144 


0.14226 


-0.035l^ 


0.293+5*0 


0.334+19 


8.6/5 


0.12100 


0.14144 


0.14262 


-0.035lJ^ 


0.297+1^ 


0.3221^3 


8.8/5 


0.12100 


0.14226 


0.14144 


-0.014+^0 


0.318+5'i 


0.323^6^ 


8.6/5 


0.12100 


0.14226 


0.14226 


-0.022^^5 


0.324l|^ 


0.323+^^ 


6.2/5 


0.12100 


0.14226 


0.14262 


-0.028ii 


0.330ti 


O.SUtlr 


5.6/5 


0.12500 


0.14144 


0.14144 


-0.1181^ 


0.190+^ 


0.353+6^ 


16.4/5 


0.12500 


0.14226 


0.14144 


-0.1041? 


0.214t^ 


0.333t^^ 


6.6/5 


0.12900 


0.14144 


0.14144 


-0.188^4 


o.iosti 


0.365t^^ 


12.6/5 


0.12900 


0.14144 


0.14226 


-0.190l^ 


0.109+^ 


0.361^8° 


7.5/5 


0.12900 


0.14144 


0.14262 


-0.190l^ 


0.112+^ 


0.353ti2 


10.8/5 


0.12900 


0.14226 


0.14144 


-0.177l^ 


0.126+? 


0.342l^^ 


6.2/5 


0.12900 


0.14226 


0.14226 


-0.182+^° 


0.129+?o 


0.339^9^ 


6.9/5 


0.12900 


0.14226 


0.14262 


-0.18611^ 


0.133+15 


0.335t?6 


12.4/5 


0.13300 


0.14144 


0.14144 


-0.240l^ 


0.045+3 


0.374+^^ 


8.5/5 


0.13300 


0.14226 


0.14144 


-0.233l^ 


0.057l^ 


0.35lt^^ 


7.5/5 



TABLE IV. Results of constant fits to T2{q , mp; my). 
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ri(m 


q) = a + bruq 




Ti(mg) = c 




Hh 


Kg 


a 


b 


xVd.o.f. 


c 


xVd.o.f. 


0.121 


0.14144 


0.299ttl 


-0M2ttll 


0.04/1 


0.2811^^ 


0.3/2 


0.121 


0.14226 


0.2891^^ 


-0.326tJli^ 


0.05/1 


0.27ll?^ 


0.1/2 


0.121 


0.1419 


0.29311^ 






0.275+1^ 




0.129 


0.14144 


0.336+^^ 


-0.058ti? 


0.1/1 


0.333+^^ 


0.1/2 


0.129 


0.14226 


0.330li? 


-0.288t|?i 


0.2/1 


0.316^9^ 


0.8/2 


0.129 


0.1419 


0.333+1^ 






0.324+^^ 





TABLE V. Extrapolation of Ti{q^ = 0), from a single pole fit, to the chiral limit, where 
Ti is assumed either to have a linear dependence on the pole mass of the light quark, or to be 
independent of the pole mass. Kstmnge = 0.1419 corresponds to the physical strange quark mass 
from determining the mass of the K on this lattice. 





T2(g^=0; mq) = a + brUq 




T2((?2=0;mg) = c 




I^h 


K,g 


a 


b 


xVd.o.f. 


c 


xVd.o.f. 


0.12900 


0.14144 


0.363+il 


-0.31lt^^o 


0.1/1 


0.348t^^ 


0.6/2 


0.12900 


0.14226 


0.337+3? 


-0.190t^^^ 


0.1/1 


0.329+9^ 


0.3/2 


0.12900 


0.14190 


0.349+^? 






0.337+^^ 




0.12100 


0.14144 


0.323t^3 


-0.406ti^ 


0.01/1 


0.304tiJ 


0.3/2 


0.12100 


0.14226 


0.311+4? 


-0.531™ 


0.1/1 


0.289+1^ 


0.6/2 


0.12100 


0.14190 


0.317+i 






0.296+1^ 





TABLE VI. Extrapolation of r2(g^=0), from a single pole fit, to the chiral limit, where T2 
is assumed either to have a linear dependence on the pole mass of the light quark, or to be 
independent of the pole mass, ^strange = 0.1419 corresponds to the physical strange quark mass 
from determining the mass of the K on this lattice. 
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T2iq?naxi^q) = a + bnig 




T2iq'iax-,mg) = c 




K-h 


Kg 


a 


b 


xVd.o.f. 


c 


xVd.o.f. 


0.12100 


0.14144 


0.331+?9 


0.492+348 


0.4/1 


0.353+^^ 


1.9/2 


0.12100 


0.14226 


0.327t?7 


-0.026+378 


1.9/1 


0.325+^^ 


1.9/2 


0.12100 


0.14190 


0.328+1^ 






0.337+g^ 




0.12900 


0.14144 


0.370t^;^ 


-0.136^1^8 


0.9/1 


0.3631^3 


1.1/2 


0.12900 


0.14226 


O.349t?0 


-0.136tlll 


0.7/1 


0.341+^3 


0.8/2 


0.12900 


0.14190 


0.358+i!^ 






0.351+g3 





TABLE VII. Extrapolation of T2{qmax) to the chiral limit, where T2 is assumed either to have 
a linear dependence on the pole mass of the light quark, or to be independent of the pole mass. 
^'Strange = 0.1419 Corresponds to the physical strange quark mass from determining the mass of the 
K on this lattice. 



I^h 


mK*/mp 


T2{qlax) 


2 / 2 
qmax/''^Psi 


^2(0) (a) 


Ti(0) (b) 


Ti(0) (c) 


?2(0) (d) 


T2(0) (c) 


0.13300 


0.59l^ 


0.362+^^ 


0.08+} 


0.333+^^ 


0.35611^ 


0.3531^9 


0.359l^^ 


0.36ll^4 


0.12900 


0.48+2 


0.353+6^ 


0.15+1 


0.301+7^ 


0.3241^9 


0.326^7^ 


0.339+9^ 


0.352+f 


0.12500 


0.421^ 


0.346+^^ 


0.20+1 


0.276l7^ 


0.29811? 


0.298+1^ 


0.318^1^ 


0.342l^^ 


0.12100 


0.37li 


0.339+7^ 


0.26li 


0.252+^^ 


0.278+14 


0.276li3 


0.298li^ 


0.332+f 




0.1692+i 


0.269+^^ 


0.51+^ 


0.112l^ 


0.124+20 


0.159+34 







TABLE VIII. Comparison of results from different methods of extracting Ti^2('7^=0). The 
last row indicates the final extrapolation to the physical regime niK'/'mB-, with results of the 
three methods of extracting Ti^2{q'^= 0','m'B',rnK*) at the feet of columns labelled (a), (b) and (c). 
Column labels: (a) pole form, with 1+ mass from two-point functions [estimated from Bgi mass 
for final extrapolation], (b) dipole form, with rrip scaling for final extrapolation, (c) pole form, 

1/2 

with rrip scaling for final extrapolation, (d) pole form, with mass determined from pole fit, (e) 
constant form factor. 
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FIGURES 




FIG. 1. An example of a penguin diagram contributing to the decay b -^ sj. 
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FIG. 2. A typical plot of Ti(g^=0; mp; my) vs. time. From the application of the time reversal 
operator, it can be shown that only the real component of Ti is non-zero. 
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FIG. 3. Im(r2), for a typical momentum used. From the application of the time reversal 
operator, it can be shown that only the imaginary component of T2 is non-zero. 
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FIG. 4. A typical plot of T2(g^Q^; mp;r7T,y) vs. time. 
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FIG. 5. Ti{q'^), using a pole fit. The dotted lines represent the 68% confidence levels of the fit 
at each q'^. 
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FIG. 6. Chiral extrapolation of Ti(g^=0). The dotted lines indicate the 
of the fit. arrig ng^f is the lattice pole mass. 
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FIG. 7. T2(g^), with a pole fit. The dotted lines represent the 68% confidence levels of the fit 
at each q^ . 
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FIG. 8. Chiral extrapolation of T2(g, 



max) 
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FIG. 9. The ratio T1/T2 at q^=0 for dipole/pole and pole/constant fits. 
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FIG. 10. The ratio Ti{q =0;mp;mK*)/T2 {q =0;mp;mK*) with Ti fitted to a dipole form 
and TJ" ^ extrapolated from T2{q'^n^) using a fitted 1^ mass. 
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FIG. 11. Extrapolation of T2(g4ax) ^^ ms, assuming HQET. The quantity plotted is T2, defined 
in the text, which is equal to T2{qmax) ^^ "^-B- 
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FIG. 12. Ti extrapolation, for N = 1/2,3/2 (Points at rriK'/'mB displaced slightly for clarity). 
Ti is defined in the text and agrees with Ti(g^=0) at niB for both N = 1/2, 3/2. The CLEO point 
is obtained from the CLEO measurement of BR{B -^ K*^) as explained in the text. 
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